Three different routes from the directed Ising to the directed percolation class 
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Scaling nature of absorbing critical phenomena is well understood for the directed percolation 
(DP) and the directed Ising (DI) systems. However, a full analysis of the crossover behavior is still 
lacking, which is of our interest in this study. There are three different routes from the DI to the DP 
classes by introducing a symmetry breaking field (SB), breaking a modulo 2 conservation (CB), or 
making channels connecting two equivalent absorbing states (CC). Each route can be characterized 
by a crossover exponent, which is found numerically as <p = 2.1 ± 0.1 (SB), 4.6 ± 0.2 (CB), and 
2.9 ± 0.1 (CC), respectively. The difference between the SB and CB crossover can be understood 
easily in the domain wall language, while the CC crossover involves an additional critical singularity 
in the auxiliary field density with the memory effect to identify itself independent. 

PACS numbers: 64.60.Ht,05.70.Ln,89.75.Da 



I. INTRODUCTION 

Critical behavior of absorbing phase transitions can be 
categorized by universality classes [l|. The most well- 
known example is the directed percolation (DP) class 
which may be (loosely) defined by a classification scheme 
known as the "DP conjecture" 0, S S]- The conjec- 
ture states that models exhibiting a continuous transi- 
tion into a single absorbing state should belong to the DP 
class unless symmetry or conservation is involved. The 
DP conjecture is not complete in two folds. First, some 
models with infinitely many absorbing states such as the 
pair contact process @ also belong to the DP class, at 
least in its stationary property, unless there is no clear- 
cut symmetry among absorbing states [a, 0, la]- Second, 
the conjecture itself does not resolve the controversy re- 
lated to the pair contact process with diffusion (PCPD) 

However, this conjecture is useful enough to trigger 
a search for a non-DP universality class. One can find 
several universality classes in Ref. [l[ which do not meet 
the constraint of the DP conjecture. Among them, this 
paper interests in the system with two equivalent absorb- 
ing states (or in general two equivalent groups of many 
absorbing states) [H, [H, [13, [3, El. The Ising- like Z2 
symmetry between two absorbing states which is some- 
times manifest as a modulo 2 conservation of domain 
walls in one dimension renders the system to exhibit a 
non-DP critical behavior. This universality class has sev- 
eral names such as the parity conserving [la . [l6[ , the di- 
rected Ising (DI) [13, [20, [2II which will be used to refer 
to this class in this paper, the generalized voter class [23 
and so on. In higher dimensions, all these classes become 
distinct. 

The Ising symmetry (or the modulo 2 conservation) per 
se is not enough to force the DI critical behavior. It re- 
quires an infinite dynamic barrier between two equivalent 
(group of) absorbing states [0, [20] , which is analogous to 



the free energy barrier between two ordered states in the 
equilibrium Ising model. A state near one absorbing state 
cannot evolve into the other absorbing state within a fi- 
nite number of successive local changes. In other words, 
a domain wall (frustration) in a configuration generated 
by pasting two absorbing states cannot disappear by it- 
self, so the system never becomes absorbing with a single 
frustration [6l.[l(j|. 

As an example without this additional feature, con- 
sider a one-dimensional system of diffusing particles with 
pair annihilation (2A -^ 0) and branching of two parti- 
cles by a triplet {3 A -^ 5 A). We may map this particle 
model onto the ferromagnetic Ising spin model by inter- 
preting a particle as a domain wall of the Ising system as 
follows: 



AID^U 

AA^m 

0AAA0 ^ AAAAA 



TiT^TTT, 
iiTiTT^iTitiT, 



(1) 



where A (0) stands for a particle (vacancy) and t (i) rep- 
resents a up (down) spin. Clearly the number of domain 
walls is conserved modulo 2 and there are two absorbing 
states in the spin system (all spins are up or down) which 
arc equivalent. However, the states with only one (dif- 
fusing) particle are also absorbing in the particle model. 
In the spin language, these extra absorbing states can be 
obtained by connecting two equivalent absorbing states 
(for example, . . . TTT-Li-L ■ • Oi which violate the infinite 
dynamic barrier requirement for the DI universality class 
mentioned above [y, [20| . It is numerically shown that this 
system shares the critical behavior with the DP not the 
DI [lOl . [23| . while a similar system with 2 A -^ and 
2 A -^ iA belongs to the PCPD class [H]. It is claimed 
that a long-term memory plays a crucial role in differen- 
tiating these two universality classes |ll| . 

According to the above criterion for the DI class, 
one can conceive three possible routes along which the 



DI critical behavior crosses over to the DP in one di- 
mension. Introducing a symmetry breaking field (SB) 
[23, S S [13, [23 , allowing dynamics which breaks the 
modulo 2 conservation (CB) [2^, [30[, or making channels 
connecting two equivalent absorbing states (CC) [y, Q 
drives the system into the DP class. But these three 
routes bring about different characteristics in the absorb- 
ing states. The SB route maintains the absorbing states 
but breaks the probabilistic balance between them. The 
CB route (for example, by adding A — > or 2A) in- 
validates the mapping between the spin and the parti- 
cle model, and ends up with only one absorbing state 
(vacuum). The CC route is accompanied by additional 
(infinitely many) absorbing states connecting two equiv- 
alent absorbing states with both the symmetry and the 
mapping sustained. 

To quantify the difference of these routes, we measure 
the crossover exponents governing the crossover behav- 
ior. In Sec. |TT1 one-dimensional models are introduced 
belonging to the DI class. In Sec. IIIIl numerical results 
for the crossover exponents along three different routes 
are reported. We discuss and summarize our results in 
SecHYl 



II. DI MODELS 

We choose two DI models for studying crossover from 
the DI class to the DP class: branching annihilating 
walks with two offspring (BAW2) [H, [l^ and the in- 
teracting monomer model (IM) which is a simplified ver- 
sion of the interacting monomer-monomer model [8[ . We 
study the macroscopic absorbing phase transitions [23[ 
inherent in the model. The suitable order parameters 
will be defined when we explain each model in detail. For 
convenience, the order parameter for different models will 
be denoted by p throughout the paper. We only study 
one dimensional lattice systems with periodic boundary 
conditions. 

For the simulation of the BAW2, we take the dynamic 
branching rule of Ref . [29| . The simulation algorithm is 
detailed as follows: At first, we choose a particle resid- 
ing at a lattice point randomly among, say, Nt particles 
present at time t. It may hop to a randomly chosen near- 
est neighbor site with probability p. If the target site is 
already occupied, two particles are removed immediately. 
With probability (1 — p), this selected particle branches 
two particles along one of two directions. If an offspring 
is created on a site already occupied, both particles anni- 
hilate immediately. After the above steps, time increases 
by l/Nf. Since the absorbing state is vacuum, the order 
parameter p is the particle density. At criticality, the 
density is expected to decay as p{t) ~ t^^ with the DI 
critical exponent S ~ 0.285 [13 • Figure [1] locates the 
critical point at Pc = 0.510 35(5) which is more accurate 
than that obtained in Ref. [29J . 

The IM is a lattice model with two possible states at 
each lattice point; either vacant or occupied by a particle. 
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FIG. 1: (Color online) Semilogarithmic plot of p{t)t vs t 
near criticality for the BAW2. The critical point is found 
to be Pc — 0.510 35(5) with the figure in the parentheses as 
uncertainty of the last digit. 



The dynamics begins with the adsorption attempt at a 
randomly selected vacant site. When both of nearest 
neighbor (nn) sites are vacant, a monomer is adsorbed 
with rate 1 . In case only one of the nn sites is occupied, 
the monomer is adsorbed on the substrate with rate A, 
then reacts and leaves the substrate with the neighboring 
particle (^-l-A — s- 0) instantaneously. When both nn sites 
are already occupied, the adsorption attempt is rejected. 
The dynamics on the substrate can be summarized as 



VVV -U VAV, 

^yn ^ VVV 

VVAJ -^^^^^ 



(2a) 
(2b) 



where A (V) means the monomer occupied (vacant) site. 

The absorbing states of the IM are characterized by the 
"antiferromagnetically ordered" state (. . . AVAVAV . . .) 
and have the Z2 symmetry (poisoning even or odd sites 
by monomers). Starting with initial configurations with- 
out any AA pair, the proper order parameter p is the 
density of VV pairs. 

The mapping of the dynamics of the IM to that of the 
contact process with the modulo 2 conservation [CP(2)] 
[3l| is possible, though there are some subtleties which 
make these two models a little bit different. If we assign 
a domain wall excitation, say X, in the middle of the 
VV pair, the domain wall dynamics is exactly the same 
as the CP(2) rules. A pair annihilation in the CP(2) cor- 
responds to an adsorption process of a monomer as in 
Eq. (Pa|) and a pair branching in the CP(2) to an anni- 
hilation process Eq. (j2bp . Note that the reaction (j2bp 
always generates two domain walls (two VV pairs), be- 
cause an AA pair is not allowed in the IM. However, there 
is no one-to-one correspondence between configurations 
of the CP(2) and the IM, in general. For example, a 
domain wall configuration A0A is not possible in the 
IM by definition, but possible in the CP(2) in principle. 
Nevertheless, if we start with the fully occupied domain 
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FIG. 2: (Color online) Semilogarithmic plot of p{t)t vs t 
near criticality for the IM. The critical point is found to be 
Ac = 0.8930(1). 



wall state, such a configuration cannot be generated by 
dynamics and thus the map ping becomes exact. In this 
case, the parameter q in Ref. [31[ for the CP(2) is related 
to 1/(1 + 2A) of the IM with suitable time rescaling. This 
mapping guarantees that the IM should belong to the DI 
class which is confirmed by Fig. [2 We found the critical 
point of the IM to be Ac = 0.8930(1) and correspond- 
ingly the more accurate critical point of the CP(2) as 
q^ = 0.358 94(3) than that in Ref. [3l|. 



III. THREE ROUTES TO THE DP CLASS 

To study the crossover from the DI class to the DP 
class, we introduce dynamics with rate w which breaks 
either the symmetry or the mod-2 conservation. The 
crossover scaling ansatz of the order parameter at criti- 
cality is [ii,[3l] 



p(i;A,'u;) =t-''F{ll''H,w^'H), 



(3) 



where A is the deviation of a tuning parameter from the 
DI critical point. The meaning of A and w will be made 
clearer in the context of each crossover model. When 
0, the scaling function in Eq. ^ should describe 

5 and v\\ are the critical 
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the DI critical behavior. Hence, 
exponents of the DI classes; 5 
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3.2, 



respectively [19|. The crossover exponent is defined as 
(j) = v\\/^^\\ which can be measured from observing how 
the phase boundary behaves near the DI critical point 



hiAc(«;) 

hm — 

■u)^o+ Inw 



(4) 



where Ac(w) is the distance from the DI critical point at 
w = to the DP critical point at nonzero w. 

First, consider the conservation-breaking (CB) route 
by introducing the mod-2 conservation breaking dynam- 
ics to the BAW2. When a particle tries to branch off- 



spring, the number of offspring is 2 (1) with the proba- 
bility 1 — w (w) . This crossover model with finite w will 
be referred to as the BAWCB. 

Second, the symmetry-breaking (SB) route is con- 
structed by introducing the symmetry breaking dynam- 
ics to the IM in such a way that the adsorption rate 
Eq. ((2a|) at even-indexed sites is different from that of 
odd-indexed sites. To be specific, we modify adsorption 
dynamics without affecting the dcsorption events as 



VVoV 



VAoV, VV,V 



VA,V, W) 



where the subscripts o and e refer to an odd- and even- 
indexed sites, respectively. When < u; < 1, this sys- 
tem prefers the absorbing state of poisoning monomers at 
even-indexed sites, which breaks the Z2 symmetry of the 
IM. We will refer to the model with dynamics Eqs. ((2k^|) 
and jSEI as the IMSB. 

Finally, the channel-connecting (CC) route can be con- 
structed by connecting two equivalent symmetric absorb- 
ing states through infinitely many absorbing configura- 
tions Q . In the context of the IM model, this route can 
be studied by allowing an adsorption without a pair re- 
action probabilistically: 



AW 



AAV, VVA 



VAA. 



m) 



The model with dynamics of Eqs. ^^, (l2b|, and ((20) 
will be called as the IMCC. As an adsorption attempt 
at a vacant site between two monomers (AVA) is still 
rejected, any configuration without a VV pair is absorb- 
ing and the order parameter is still the VV pair den- 
sity. In contrast to the IM, the IMCC has infinitely 
many absorbing states characterized by the AA pair den- 
sity (auxiliary field density). By pasting the two anti- 
fcrromagnetically ordered absorbing states, we find ei- 
ther . . . VAVAAVAV . . . which is already absorbing, or 
. . . VAVAVVAVAV . . . which can evolve back to an ab- 
sorbing state by adsorbing a monomer at cither site of 



TABLE I: Critical point values Ac(u>) of the IMSB and IMCC 
and Pc{w) of the BAWCB for various w's. The numbers in 
the parentheses indicate the uncertainty of the last digits. 

w Xcjw-IMSB) \c{w-lMCC) pc(w; BAWCB) 

0.8930(1) 

10"^ 

2 X 10"^ 

5 X 10~^ 

10-^ 

2 X lO"'* 

3 X 10"^ 
5 X 10-^ 

10"^ 
2 X 10~^ 

4 X 10"^ 



0.9125(2) 
0.9199(1) 
0.9303(2) 
0.944 75(5) 



0.8930(1) 
0.9119(1) 

0.92545(5) 
0.9340(1) 
0.944 95(5) 

0.965 35(5) 
0.9842(1) 



0.510 35(5) 
0.4931(1) 
0.490 30(5) 
0.4858(1) 
0.481 65(5) 
0.4767(1) 
0.473 45(5) 
0.468 90(5) 
0.461 75(5) 
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FIG. 3: (Color online) Log-log plot of Ac(it)) vs w for three 
different crossover models. For graphical clarity, we shifted y 
axis with arbitrary scale. Each line whose slope is l/(j> shows 
the fitting result for the corresponding phase boundary. 
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FIG. 4: (Color online) Long-time behavior of the auxiliary 
field density pAA for small w's. The exponent fop — 0.160 
is the temporal decay exponent of the order parameter in the 
DP class. Simple extrapolation lines are given by dotted lines. 



the VV pair. Hence there is no infinite dynamic barrier 
between absorbing states and the system leaves the DI 
class into the DP class. 

Now we are equipped with three different models which 
can show the crossover behavior from the DI class to 
the DP class. As explained above, the crossover ex- 
ponent (f> can be deduced from the DP phase bound- 
ary near the DI critical point; see Eq. ([4]). For each 
model, the critical points at finite w are located numer- 
ically using the DP values of the critical exponents (see 
Table [J). The critical points for the IMSB and IMCC 
are denoted by \c{w) and for the BAWCB by pdw), re- 
spectively. To measure the crossover exponent, we define 
Ac{w) = \\ciw) - Ac(0)|/Ac(0) for the IMSB and the 
IMCC, and similarly Ac{w) = \pc{w) ~ Pc{0)\/pc{0) for 
the BAWCB. 

In Fig. [21 the crossover exponents are estimated for all 
three different routes. Along the SB route, we estimate 
(j) = 2.1(1) which is consistent with previous estimates: 
= 2.1(1) d^, 2.24(10) m, and recently 1.9(1) [H. 
Along the CB route, our estimate is 6 = 4.6(2) which is 
close to the recent result (f> = 4.8(2) [33|. Along the CC 
route, we find (f) = 2.9(1), which is clearly different from 
those for the SB and CB routes. 



IV. 



DISCUSSION AND SUMMARY 



The crossover behavior from one fixed point to another 
docs not reflect the properties of both fixed points, in 
general. Rather, strictly speaking, it is related to one 
fixed point and its crossover operator which forces the 
system to crossover to the other fixed point [ij, [S^, [3J| . 
In many cases, the crossover route is unique between two 
fixed points, so is the crossover operator, which leads to 
one unique crossover exponent between two universality 
classes. 



However, in the case studied here, we found that there 
are three different routes for the crossover from the DI 
to the DP class, which yield all three different crossover 
exponents. It indicates that the crossover operator as- 
sociated with each route should be different from each 
other. The difference between the SB and CB routes has 
been noticed earlier [20| . Summarizing it in the domain- 
wall (particle) representation, the crossover operator in 
the CB route is a single particle annihilation or creation 
operator in the action, while that in the SB route is a 
nonlocal string operator (global product of particle num- 
ber operators) [20| . Hence it is not surprising to find that 
the two crossover exponents are different. 

The crossover behavior along the CC route is tricky. 
First, note that the Z2 symmetry between the odd- 
indexed and even-indexed sites is not broken in this route, 
which implies that the CC route is different from the SB 
one. Second, in terms of domain walls, the CC process, 
Eq. ([^ . clearly breaks the mod-2 conservation in the 
number of VV pairs (primary field, say, X), because of 
which one may naively guess the CB-type crossover along 
the CC route. However, unlike the ordinary CB models, 
the CC process generates an infinite number of absorb- 
ing configurations with AA pairs (auxiliary field, say, Y) 
which give the feedback effect to the domain wall dynam- 
ics. Consider the dynamics of the IMCC in terms of two 
fields X and Y . The reactions, Eqs. (Pa|) and (|2b[) . can 
be rewritten as 2A ^ and X -^ iX, while Eq. ((20) 
is equivalent to X —^ Y . The presence of AA pairs 
makes it possible to allow another reaction of XY -^ 2X 
(AVVAA ^ AVVVA) in combination with Eq. ((2b)) . 
The last two reactions involving the auxiliary field repre- 
sent the feedback which generates memory effects on the 
primary field. 

The emergence of the auxiliary field density and the 
feedback mechanism may be responsible for the new 
crossover behavior along the CC route. To see how the 
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FIG. 5: (Color online) Log-log plot of pAA in the steady state 
versus w. The slope of the dotted line is 0.645. 



auxiliary field may affect the crossover, we measure the 
auxiliary field density (natural density) paa{w) along the 
DP critical line. We expect its long-time temporal behav- 
ior as pAAiw^t) =paa{w) — ht~^°^ starting from the A 
particle vacuum [8|, |35| , where (5dp is the temporal de- 
cay exponent of the order parameter in the DP class. 
Figure [H confirms our expectation and the asymptotic 
values of paa is plotted against w in Fig. [5l It is clear 
that paa{^) = at the DI critical point. Near w — Q, 
the power-law singularity is found as paa ~ w" with 
a = 0.65(2). We test the universality of the CC crossover 
scaling as well as the singularity of the auxiliary field 
density by investigating various different models like a 
CC variant of the interacting monomer-dimer model [36[ 
and a two-species particle model defined as below. It 
turns out that these models form one crossover univer- 
sality class with the same CC-type crossover exponent (p 
and the auxiliary field density exponent a [36|. 

To examine the exclusive role of the feedback (mem- 
ory) effect, we introduce a two-species (X, Y) particle 
model where the feedback process can be directly con- 
trolled. The primary particles X can diffuse, while the 
secondary particles Y are immobile. Each species parti- 
cles are hard core particles, but a simultaneous occupa- 
tion by different species at a site is allowed. Therefore 
each site can be in one of four possible states such as 0, 
X, y, and XY. 

The reactions dynamics is symbolically summarized as 
2X ^^, X ^iX, X -^ y, and XY ^ Y, 2X, or 3X. 
To be more specific, the evolution rule is given as follows: 
First, choose one of X particles randomly. With proba- 
bility w, the chosen X spontaneously mutates or annihi- 
lates as X ^ y or XY -^ Y. With probability p — w, 
the X hops to one of its neighboring sites. Whenever two 
X particles attempt to occupy the same site, they anni- 
hilate immediately like in the BAW2. With probability 
1 — p, the X can branch particles in the neighborhood 
as X ^i- X + 2X with the dynamic branching rule [2^ . 
To control the feedback effect, we introduce a parameter 
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^ feedback r = 0.5 

o no feedback 
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FIG. 6; (Color online) pviw) vs w along the phase boundary 
for the two-species model with (r = 1 and 0.5) or without 
(r = 0) feedback processes as well as the variant of the contact 
process (see text) in double logarithmic scales. The slopes of 
the straight lines are 0.70, 0.71, 0.80, and 1 from above. 



r in the branching process when a Y resides along with 
the X at the same site, such that XY -^ X + 2X with 
the relative rate (1 — r) and XY -^ X + X with r. 

At w = 0, there is no spontaneous annihilation (or 
mutation) process and the model becomes exactly the 
same as the BAW2 after some transient period needed 
for removing all Y particles via branching processes. At 
finite if, we expect the DP scaling with a finite density 
of y particles; py{w) > because X particles gener- 
ate y particles via mutation processes, which cannot be 
fully ehminated via branching processes. So the CC-type 
crossover is expected in general. 

However the r = point is special. There, the pres- 
ence of y particles never affects the dynamics of X par- 
ticles, so there is no feedback mechanism to alter the 
primary particle density through the secondary parti- 
cles. The primary field dynamics is basically identical 
to the BAWCB model introduced in Sec. HIT] and the CB- 
type crossover scaling should appear at r = along with 
nonzero py{w). At finite r, the X particle dynamics is 
affected by the presence of Y particles, which generates 
the memory effects on the X particle density. Hence r 
can be regarded as the feedback controlling parameter 
and w as the crossover parameter. 

We performed numerical simulations for various values 
of r. For nonzero r, we found that the secondary particle 
density behaves as py{w) ~ w" with a — 0.70(3) and the 
crossover exponent cj) — 3.1(2), both of which are consis- 
tent with the CC crossover results for the IMCC model 
within numerical errors. In contrast, at r = the CB 
crossover is found as expected along with a ~ 0.80 dif- 
ferent from the CC crossover case; see Fig. |6l This leads 
to the conclusion that the feedback process (memory) is 
crucial in establishing the CC-type crossover universal- 
ity class and affects the singularity characteristic of the 



auxiliary field density. 

It is interesting to note that the numerical values of a 
suggest a simple conjecture that a = 1 — 1/0 for the DI to 
the DP crossover with an emerging auxiliary field. The 
mean field analysis leads to a = and cj) = 1 (see Ap- 
pendix A) which is also consistent with this conjecture. 
The nontrivial singularity in py is related to the nontriv- 
ial crossover from the DI to the DP. Then, it is natural to 
ask what will happen if there is no nontrivial crossover. 
To answer this question, we study a two-species version 
of the contact process (CP) with the secondary parti- 
cles; X ^ 9, X ^ 2X, X ^ Y, and XY -> Y, X, 
or 2X with the crossover parameter w controlling the 
X ^ Y process and the feedback parameter r control- 
ling the XY — > X. This model belongs to the DP class 
for any value of w and r, including w = 0. Therefore 
there is no "true" crossover in this model. We found the 
trivial linear phase boundary (Ac(w) '^ w) as expected 
and also the trivial behavior of pY with a = \ regard- 
less of the presence of the feedback processes (see Fig. [6]) . 
These results are consistent with the mean field results 
described in the Appendix A. 

Up to now, we have been only interested in the emer- 
gence of the auxiliary field for a finite crossover parameter 
w which breaks the mod(2) conservation for the primary 
particles. One may consider more general cases where the 
auxiliary field density py is finite even at the DI critical 
point [w = 0). First, we studied the crossover behav- 
ior to the DP models with a continuous variation of py 
as w increases. We found that py{'w) — py(0) ~ w°-^^ 
and the crossover scaling belongs to the CB class |36| . 
This implies that the vanishing py (not just singular- 
ity near w = 0) is another important ingredient in the 
CC crossover. Note that py{w) seems to behave in the 
same way as the phase boundary Ac(u') ^ w^''^ with 
l/(f) ~ 0.22 for the CB crossover. 

Second, there may be a discontinuous drop in py in 
the crossover to the DP with a single absorbing state 
{py = 0). This case may be compared to the crossover 
from the DP with finite py (infinitely many absorbing 
states) to the DP with py — and similarly from the 
DI with finite py to the DI with py = [3J| . As under- 
stood in Ref. [3J], this crossover may be characterized 
by a discontinuous jump in the phase boundary at the 
DI critical point for the excitatory route or a continuous 
phase boundary for the inhibitory route. Along the in- 
hibitory route, the CB crossover is observed again (not 
shown here) and we conclude that the discontinuity in the 
auxiliaryfield density does not provide any new crossover 
scaling |36| . 

To summarize, we studied three routes of the crossover 
from the DI to the DP class; symmetry breaking (SB), 
mod-2 conservation breaking (CB), and channel connect- 
ing (CC) routes. These three routes are characterized by 
three different crossover exponents. The difference be- 
tween the SB route and the CB route is clear because 
the symmetry breaking field can be interpreted as the 
spatial non-local operator affecting the domain wall dy- 



namics, while the conservation breaking corresponds to 
the local operator. The CB route and the CC route share 
some common features, but the feedback memory effect 
(or temporal non-locality) makes the CC crossover dis- 
tinct from the CB crossover. We also found the universal 
exponent a = 0.65(2) which describes the singularity of 
the vanishing auxiliary field in the CC route. From the 
numerical study of the feedback controlling model and 
its mean field theory, we conjectured that a = 1 — l/(j). 
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APPENDIX A: MEAN FIELD THEORY 

This appendix provides the mean field theory for the 
feedback controlling model introduced in Sec. IIVI Here 
px and py are the densities of X and Y particles, respec- 
tively and pz is the density of sites where both X and 
Y particles reside. Then one may easily write down the 
mean field equations as 



Px_ 

PX 



2q — w — 2(1 + q — w)px — rq 



Pz 
PX 



- 2pz im 



py = 


= wpx - {q + w)pz 








(A2) 


PZ 
PZ 


= -l + ipy-2pz) 


[l + q- 


■.Px 

-w) 

pz 


- rq 


, (A3) 



where q = 1 — p. Since pc{w = Q) = 1 in the mean field 
theory, q is the same as A in the text. 

In the active phase, the steady state density can be 
calculated by setting px = and so on. The superscript 
s in the following indicates the steady state density. From 
Eq. (|A2p . we get in the active phase 



Px 
P% 



= 1 



q_ 

w 



(A4) 



Hence Eq. 
that 



Py 



3]) gives the steady state density of py such 



2p| 
2p| 



(1 



w)PxIp% 
w 



rq 



9(1 



+ q) +w{l — w) 



(A5) 



In the active phase, the (stable) steady state density 
Py is determined uniquely, so we define the natural den- 
sity at criticality by taking the limiting process from 
the active side as lim + Py(g). From Eqs. (jAip and 



4|), the critical point is determined by the equation 
rq + (2(7 ~ w){\ + q/w) = which reads 



g, = |{V8 + (l-r)2-(l-r)}. 



(A6) 



For any value of r(< 1), qc ~ Ac ~ w which gives 0=1. 
Since p^ = at criticality, we get py ~ w/q^ '^ w° 
from Eq. (jASp . yielding a = 0. In low dimensions, ^c 
is expected to renormalize such that qc ^ Ac ^ w^^'^, 
which leads the conjecture a = 1 — l/(f>. 

For comparison, we study the two-species version of the 
contact process with the secondary particles introduced 
in Sec. IIVI The mean field equations are 



Px = {2q-l)px -qpx -rqpz{l- px), (A7) 
PY = wpx - {q + w)pz, (A8) 

pz = -pz + q{pY - pz){px - rpz), (A9) 



where q = 1 —p. Then one may easily show that the criti- 
cal line and the critical steady-state density of Y particles 
are given as for small w 



Ac 
Py 



r 
2^' 



9c + 9c(l - r)w 



Aw, 



(AlO) 



where Ac = qc{w) — 1/2. As expected, we find the trivial 
phase boundary and also the trivial value of a = 1, which 
is consistent with the simulation results shown in Fig. [51 
Moreover, the renormalization cannot generate a singular 
behavior in the denominator of Eq. (jA10[) . so a is always 
expected to be 1 in all dimensions. 
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